Spheroidal nanoparticles as thermal near-field sensors 
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We suggest to exploit the shape-dependence of the near-field heat transfer for nanoscale thermal 
imaging. By utilizing strongly prolate or oblate nanoparticles as sensors one can assess individual 
components of the correlation tensors characterizing the thermal near field close to a nanostructured 
surface, and thus obtain directional information beyond the local density of states. Our theoretical 
considerations are backed by idealized numerical model calculations. 
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I. INTRODUCTION 



While profound theoretical understanding of thermally 
generated electromagnetic fields close to material sur- 
faces has already been obtained some time ago [H-ll], 
accurate experimental studies of such fluctuating ther- 
mal near fields at distances on the order of or even sig- 
nificantly shorter than the thermal wave length from a 
sample's surface have only recently become possible. Hu 
et al. have measured the radiative heat fiux between 
two glass plates spaced by a micron-sized gap [1], and 
reported near-field heat transfer exceeding the far-field 
limit set by Planck's law of blackbody radiation. Chen 
and co-workers have studied the near-field heat flux be- 
tween a sphere and a flat substrate 0, 0|, focusing on 
the role of surface phonon polaritons, and observed heat 
transfer coefficients even three orders of magnitude larger 
than the blackbody radiation limit. Rousseau et al. have 
made precise measurements of the radiative heat trans- 
fer between spheres with diameters of 22 and 40_um at 
distances from 30 nm to 2.5 \xm from a plate [8|, and 
obtained impressive agreement withpredictions based 
on fiuctuational electrodynamics [l|, [1]. This finding is 
conceptually important since it indicates that fluctua- 
tional electrodynamics, which is a macroscopic theory, 
can be relied on at least up to the smallest distances ac- 
cessible in this experiment, without having to account 
for effects caused by a nonlocal optical response [§, [l3| . 
Moreover, a device termed Near-Field Scanning Ther- 
mal Microscope (NSThM) has been designed by Kittel 
et al. pH - fl^ This instrument allows one to record the 
near-field heat flux at probe-sample distances even down 
to a few nanometers, but still further efforts are required 
to achieve accurate calibration. A complementary set-up 
dubbed Thermal Radiation Scanning Tunneling Micro- 
scope (TRSTM) has enabled De Wilde et al. to take 
images of thermally excited surface plasmons, and to 
demonstrate spatial coherence effects in near-fleld ther- 
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mal emission [ij]. This notable growth of experimental 
research on thermal near-field phenomena is triggered, 
on the one hand, by a wide variety of possible technolog- 
ical applications, ranging from the design of nanoscale 
heaters for use in heat- assisted magnetic recording or 
heat-assisted lithography [1| to near-field thermophoto- 
voltaics [15''-^^. On the other hand, there is a need for 
further insight into basic small-scale phenomena. For ex- 
ample, at very short distances nonlocal effects ultimately 
must come into play fo', Tol , and the use of macroscopic 
fluctuational electrodynamics as a theoretical framework 
might no longer be sufficient. Thus, thermally induced 
near-field effects constitute an emerging subject which 
requires further experimental and theoretical investiga- 
tions; of particular importance is the identification of 
experimentally observable quantities against which theo- 
retical models can be tested. 



It has already been demonstrated experimentally that 
the NSThM opens up the possibility of thermal irnag- 
ing of surface structures with nanoscale resolution (20|, 
since the surface topography leaves its imprint in the lo- 
cal density of states (LDOS). In this paper we take a 
further step in this direction: We argue that one can ob- 
tain information even beyond the mere topography if the 
NSThM sensor is appropriately shaped. This is due to 
the fact that in principle not only the LDOS, but even the 
individual diagonal components of the fluctuating fields' 
correlation tensors are measurable; these components are 
singled out when employing strongly prolate or oblate 
spheroidal particles as sensors. Therefore, when scanning 
a surface at nanometer distances with such a sensor, the 
NSThM-signal does not give a one-to-one image of the 
surface profile, but rather provides information on how 
that profile affects the directional properties of the fluc- 
tuating thermal near fleld. We first sketch the theoretical 
idea underlying this proposal in Sec.|lll and then present 
in Sec. lIIII some numerical model calculations which illus- 
trate its principal feasibility. These calculations necessar- 
ily are strongly idealized, and do not refer to an already 
existing experimental set-up; nonetheless, they indicate 
a possible direction for future developments. 
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II. THE BASIC PRINCIPLE 



Consider a substrate with temperature Ts- The elec- 
tric field E(r,i) close to its surface, generated by ther- 
mal fluctuations inside the sample, then is characterized 
in terms of its correlation tensor Wjj.{r,r' ,uj), which is 
determined by an ensemble average 
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where Ej{r,Lu) are the field's Fourier components. The 
magnetic correlation tensor Wjl{r,r',uj) is defined in an 
analogous manner. According to fluctuational electrody- 
namics P, 0j , these correlation tensors are proportional 
to the Bose-Einstein function evaluated at the tempera- 
ture of the sample, 
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Now assume that a nanoparticle with electric polarizabil- 
ity tensor q^(iij) is brought into this fluctuating near field 
at a position r, such that its distance from the surface still 
remains appreciably larger than its characteristic linear 
size, while that size should be smaller than the thermal 
wavelength at Ts . Then the field introduces a dipole mo- 
ment p(i) in the particle which in its turn interacts with 
the field, such that the power 



(3) 



+ 00 



(-ic.)^tr(aE(^).t^E(^^^^^)^ 



is dissipated in the particle, resulting in its heating. The 
right-hand side of this equation still is written in a man- 
ner which shows that it does not depend on the choice of 
the particular coordinate system used for its evaluation, 
as the trace over the product of the polarization tensor 
and the correlation tensor manifestly remains invariant 
under coordinate transformations. When adopting a co- 
ordinate system in which the tensor oF'iio) is diagonal, 
with diagonal elements a? (w), this expression ^ takes 
the more familiar form 
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It is essential to keep in mind that here the diagonal 
components VF^(r,r,w) of the sample's field correlation 
are given in the principal-axis system of the particle's 
polarizability, not in a system attached to the sample's 
geometry or orientation. If all elements (w) are equal, 
as in the case of a nanosphere, the integrand simply is 
proportional to the trace of the correlation tensor, i.e., 
to the spectral energy density 
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FIG. 1: Ratios 
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; for a spheroidal 
gold nanoparticle at ojth = 10^* s^^. The particle's volume is 
kept constant when the aspect ratio _Rb/fla is varied. 



which is the case considered usually |21| - |24| . If, how- 
ever, the elements (cj) differ significantly from each 
other, the dissipated power can be dominated by indi- 
vidual components VF^(r,r,a;). This is the principle we 
are going to explore. 

For obtaining a convenient expression for the net power 
flow Ps<->p between the nanoparticle and the substrate, 
we factorize the correlation tensor in the form 



(27r) 



■W^fe'l(r,r,c.) 



e(w,rs)i?f,(r,c.). 



(6) 



so that the quantities (r, lS) carry the same dimension 
(sm~^) as the local density of states (LDOS) [25'|. Then 
two additional mechanisms have to be accounted for: On 
the one hand, as strongly emphasized by Dedkov and 
Kyasov [23'^ and by Chapuis et al. ^3] , the magnetic part 
H(r, t) of the fluctuating near fleld gives rise to the induc- 
tion of an effective magnetic dipole nioff (i) in the particle, 
and hence to a further contribution (moff(t) •H(t)) to the 
dissipated energy. On the other hand, there also is a flow 
of energy from the particle, with temperature Tp, to the 
sample. In total, this results in the formula 
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Here we assume that the principal-axis system of the 
magnetic polarizability tensor coincides with that 

of its electric counterpart. A particularly favorable case 
occurs for spheroidal nanoparticles (i.e., for rotational 
ellipsoids with two equal semi- axes i?b and a different 
third semi-axis Ra), for which the polarizabilities Oif{(j-') 
and Qf^(w) are known analytically [2a|. In Fig. [T] we plot 

the ratios r^^^ = Im(aa^^)/Im(a|^^^) of the absorptiv- 
ities along (a) and perpendicular to (b) the axis of ro- 
tation for spheroidal Au nanoparticles at the frequency 



10 s , roughly corresponding to the dominant 
thermal frequency at T = 300 K. We employ the Drude 
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FIG. 2: Cross section, and 3D plot, of the surface structure 
employed in the numerical model calculations. The orienta- 
tion of the spheroidal nanoparticle here corresponds to the 
scan of the near-field heat transfer shown in the later Fig. [S] 
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for the permittivity, with plasma frequency ujp — 1.4 • 
10^^ and relaxation rate 7 — 3.3 • 10^'^ s^^ for gold at 
room temperature, and keep the particle's volume con- 
stant when varying the aspect ratio i?b/^a- As expected, 
one observes ^ 1 for prolate, rice grain-like parti- 
cles {Rh/Ra ^1), whereas » 1 for strongly oblate, 
pancake-like particles {Rh/R& S> 1), in which the mag- 
netic fields can easily induce eddy currents. We expect 
that the model yields reliable results at least in the in- 
terval 1/5 < R\y/Ra, < 5 of the aspect ratio [27j . 



III. NUMERICAL MODEL CALCULATIONS 

For our model calculations we take a substrate which 
is perfectly flat except for a square pad with 5 nm height, 
30 nm width, and rounded edges, as depicted in Fig. [2] 
The coordinate system is oriented such that the x-y-plane 
coincides with the base substrate plane; the x- and y- 
axes being parallel to the pad's edges. We again assume 
that the permittivity of the substrate is described by the 
Drude model with parameters appropriate for gold. 

In Fig. [3] a) - c) we show the "electric" elements 
Dfi^{r,uj) as defined by Eq. evaluated for uj — wth 
at a distance of 30 nm above the base substrate plane. 
These elements have been computed by means of a per- 
turbative approach to first order in the surface pro- 
file [2^. The sum of these three quantities gives the 
electric part of the LDOS, plotted in Fig.[3]d). Whereas 
the LDOS merely yields a blurred image of the under- 
lying surface topography [l^l, the individual diagonal 
elements D^i^{r,Ld) contain interesting directional infor- 
mation. The corresponding results for the "magnetic" 
elements D^^{r,uj) are shown in Fig. U) 
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FIG. 3: a) - c): "Electric" elements Dfj;(r,aj) (fe = x,y,z) 
for the model surface depicted in Fig. (2] computed for uuth ~ 
10^^ s~^ at a distance of 30 nm above the base x-y-plane. The 
dielectric properties of the substrate are given by the Drude 
permittivity with parameters for gold, d): Electric LDOS, as 
corresponding to the sum of a) - c) . 
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FIG. 4: As Fig. [3l for the "magnetic" elements Dfi^{r,uj) 
Observe that the scales differ from those in Fig. [S] 



For understanding the message contained in these fig- 
ures, it is helpful to also consider the normalized eigen- 
vectors f*^-'^(r,u;) of the 3 x 3-matrix D^{r,u}), and the 
corresponding eigenvalues (i^-'^(r, w), together with their 
magnetic counterparts [29| . Intuitively speaking, the 
eigenvectors give the principal directions of the local field 
fluctuations; of course, these directions are determined 
by the surface profile. The eigenvalues then specify the 
strength of the partial fiuctuation associated with the re- 
spective principal direction. In compact matrix notation, 
one now has 



(9) 



where the jth column of the matrix S is given by the 
vector f''')(r,a;), the matrix is the adjoint of 5", and d 
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is the diagonal matrix with S^"^ (r, w) at the jth position 
on the diagonal line. Written out in components, this 
implies 
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The square \fjf\r,uj)\'^ obviously specifies the weight 
with which the eigenvalue d'^-* (r, ui) contributes to the 
fluctuations along the fc-direction. Thus, the individual 
diagonal elements Dff_{r,uj) express the distribution of 
the entire fluctuation strength, which in this language 
corresponds to the sum d^-'\^,^) and thus is pro- 

portional to the electric energy density ([5]), over the indi- 
vidual directions. Analogous remarks apply to the mag- 
netic fluctuations. This kind of directional information 
is, in principle, contained in the thermal near field above 
a nanostructured surface, and allows one to deduce prop- 
erties of the structure itself. The experimental task is to 
actually extract such information from given laboratory 
samples, and the theoretical task is to devise methods for 
doing so. 

To this aim, we suggest a further development of the 
NSThM: The tip of this instrument is functionalized to 
act as a thermocouple, so that the small temperature 
change resulting from the energy deposited in or drawn 
out of the tip translates into a measurable voltage [iTI - 
[isj . If the active volume could be given the shape of 
a spheroid, the shape-dependence of the heat transfer 
expressed by Fig. [T] could be exploited for assessing indi- 
vidual elements D^/^(r,w) by means of Eq. ([T]). Addi- 
tional experimental handles would then be provided by 
the choice of both the spheroid's orientation and its ma- 
terial. 

In order to demonstrate the feasibility of this con- 
cept, we compute the near-field heat transfer between 
our model gold substrate with temperature Tg ~ 300 K 
and spheroidal nanoparticles cooled down to Tp 100 K. 
For metallic nanoparticles the magnetic contribution to 
the total heat transfer generally is much larger than 
the electric one [1^ HJ] , so that one has to select non- 
metallic sensors when aiming at the electric components 
Dfi, (r , w) . Here we consider a rice grain-like particle with 
the permittivity of gallium nitride (GaN), described by 
the formula [30] 
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with parameters e^o = 5.4, S = 5.1, = 1.0 x 10^'' s ^, 
Wn = 9.1 X 10"s-i, r = 2.7 X 10^2 s-i, and 7 = 
7.1 X 10^4 s-^ the values of the inverse relaxation times 
r and 7 have been properly adjusted to the particle's 
temperature. We choose the half-axes i?a = 40 nm and 
i?b = 10 nm, so that ~ 23.4 at Wth- These parameters 
clearly fall outside the regime where the dipole model ([7]) 
can give quantitatively accurate results. We employ this 
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FIG. 5: Total heat transfer between a cooled spheroidal GaN 
nanoparticle (Tp = 100 K) with i?a = 40 nm, Rh = 10 nm, 
and axis of rotation parallel to the a;-axis, and the gold struc- 
ture sketched in Fig.[2](rs = 300 K), assuming that the parti- 
cle is scanning the structure with a constant height of 30 nm. 



model nonetheless, as computational tools for calculat- 
ing the near-field heat transfer between nanoparticles and 
structured surfaces are not yet available, and stress the 
qualitative nature of the following results. 

Figure [5] shows the total heat transfer according to the 
model (l7|), with the particle oriented parallel to the x- 
axis at a constant height of 30 nm above the base plane 
while scanning the substrate, as previously sketched in 
Fig. [21 In this case the heat flux mainly picks up the con- 
tribution proportional to Z?^2.(r, Wth)- As a consequence, 
the signal obtained here closely resembles the previous 
Fig. [3|a). The slight differences between the two figures 
visible above the pad's center are explained by the fact 
that D^^[r,u)th) contributes significantly here. 

For assessing the magnetic contributions, we select a 
pancake-like Au nanospheroid with i?a — 10 nm and 
i?b = 40 nm, giving « 8.5, and orient it along the 
z-axis. The total heat transfer then is dominated by 
D^^ (r, Wth) , as witnessed by the good agreement of Fig. [Sj 
with[4|c). Thus, here the recorded heat flux provides in- 
formation on a component of the magnetic correlation 
tensor. 



IV. CONCLUSION 

In this paper we have outlined that individual diago- 
nal components of the correlation tensors of fluctuating 
thermal fields close to the surface of a given material 
can be measured by near-field thermal imaging, if the 
sensor employed for recording the heat flux is appropri- 
ately shaped. Evidently this proposal still is somewhat 
sketchy and idealized: On the theoretical side, there is an 
urgent need for accurate computational tools for quanti- 
fying the near-fleld heat flux between arbitrarily shaped 
nanoparticles and structured surfaces; from the experi- 
mental viewpoint, the development of the required sen- 
sors poses a serious challenge. Thus, there may still be a 
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FIG. 6: Total heat transfer between a cooled spheroidal Au 
nanoparticle (Tp = 100 K) with iia ~ 10 nm, Rh = 40 nm, 
and axis of rotation parallel to the z-axis, and the gold struc- 
ture sketched in Fig.[2](rs = 300 K), assuming that the parti- 
cle is scanning the structure with a constant height of 30 nm. 

long way from our present matter-of-principle considera- 
tions to laboratory applications. 

Yet, there are at least two reasons for exploring this 
option. First, fiuctuational electrodynamics is a macro- 
scopic theory; it would be important to test its predic- 



tions even for probe-sample distances of a few nanometers 
and below, where its validity can no longer be taken for 
granted. Possible deviations form the macroscopic theory 
would then yield important novel insight into fundamen- 
tal processes in dielectrics. For this program sensitive 
observables are required; we suggest that measuring indi- 
vidual diagonal components of the near fields' correlation 
tensors, which give directional weights to the total fluctu- 
ation strength, might yield more information than their 
traces, i.e., the local density of states. Second, leaving 
theoretical considerations aside, the practical exploita- 
tion of the principles outlined in this note for nanoscale 
thermal imaging [20j would provide novel tools for the 
diagnosis and characterization of nanostructures. 
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